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ABSTRACT 

The consequences of the shear-induced alpha effect to the long-term modulation of 
magnetic activity are examined with the help of the axisymmetric numerical dynamo 
model that includes the self-consistent description of the angular momentum balance, 
heat transport and magnetic field generation in the spherical shell. We find that the 
shear contributions to alpha effect can complicate the long-term behaviour of the large- 
scale magnetic activity and differential rotation in nonlinear dynamo. Additionally we 
consider the impact secular magnetic activity variations to the secular modulations of 
the solar luminosity and radius. 
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1 INTRODUCTION 



It seems that the long-term variations of the sunspot's ac- 
tivity and phenomena similar to the Maunder minimum are 
widely presented on the late-type stars like the Sun, see 
e.g.([25||T|. Generally, it may be usefuU to distinguish the 
different kinds of time variations of magnetic activity on the 
stars for two type. One is irregular disappearence of sunspot' 
cycle similarily to what was happened during the Maun- 
der minimum at the second half of the 17-th centure, Q. 
Another kind of variations can be attributed to the quasi- 
regular grand activity cycles. The grand activity cycles is 
natural to explain as a result of the nonlinear interactions 
between magnetic fields and the main sources of the solar 
dynamo, which are the differential rotation and convection, 
([6][29][27][19][8|. The given nonlinear mechanisms may be 
responsible for the chaotic variations of magnetic activity 
on the late type stars and on the Sun, as well. Though, 
the chaotic behaviour was demostrated only in the case of 
the highly supercritical dynamos, (|6 27 8|. Currently, some 
doubts are cast upon the possibilitiy for the supercritical 
regimes in the stellar dynamo. The main constraint for fea- 
sibility of the supercritical dynamo is due to magnetic helic- 
ity conservation. However, the problem can be overcome if 
to allow the helicity outflux of the dynamo region. Another 
contraint come into the play if to consider the distributed 
dynamo. In this case the energy of the generated large-scale 
magnetic fields (LSMF) is hardly to exceed the equipartition 
level with the energy of convective fiow. These arguments 
motivate us to look for mechanisms which could induce the 
long-term variations for such a weakly non-linear dynamo 
regimes. 

In this paper we consider an effect that can provide a 
nonlinear coupling between two main sources of the solar 
dynamo, these are the differential rotation and a effect. In 
according to ( |12| |pip05| |18[ ) the shear (differential rotation) 
changes both the structure and the amplitude of the alpha 



effect. For the Sun the most important component of alpha 
effect is a^0. The estimation of shear contribution to the a^^ 
with regards fo the typical conditions in the solar CZ shows 
that it gives only about 15% of the whole magnitude. In 
the nonlinear dynamo the LSMF are included in the angu- 
lar momentum balance inside CZ. In fact the grand activity 
cycle may be considered as a result of relaxion of the pertu- 
bations angular momentum balance in CZ emerged due to 
LSMF (via the A-quenching( 6 ; 19), or the large-scale Lorenz 
force |l3l[26l [8l)). In this case the incorporation of the shear- 
induced alpha effect can complicate the relaxion process be- 
cause it supplies the additional destabilization the sources 
of magnetic filed generation due to the time-dependence of 
the alpha effect on the differential rotation distribution in 
CZ. 

Below, we demonstrate the significance of the shear- 
induced alpha-effect for the long-term evolution of magnetic 
activity using the numerical model of the axisymmetric dy- 
namo in the convective shell. The model includes the numer- 
ical solution for the dynamo equations, angular momentum 
balance and heat transport problem with reagards for the 
physical conditions that take place in the main part of the 
solar CZ. In the next part of the paper we give the short 
description of the model (more details can be found in the 
article( [l7| l). The results of numerical simulations are dis- 
cussed at the third part of the paper. 



2 THE MODEL DESIGN 

The basic physical quantities under considerations are 
the induction vector of the axisymmetric LSMF, B = 

e^B{r, 6)+curl (^-^r^^^^j^^ , the large-scale flow - differential 

rotation, V — e^r sin 6 Q (r, 9) (we neglect the meridional 
circulation, here), the mean specific entropy, S{r,d), and 
the pressure, P(r,6). The magnetic activity and differen- 
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tial rotation give rise to pertubations of the thermodynamic 
stage of SCZ. The thermodynamic deviations from refer- 
ence model (taken on the base of solar interior model from 
Stix(28l) is described with help of the mean entropy pro- 
duction equation. In following to Pipin & Kitchatinov( 20 1 
we write it as 



pT 



dS 



dt 
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Again, for details please look at paper |T7| and references 
there. The heat flux F is a sum of convective and radiative 
fluxes F = Fc -I- Frad, where in convective flux we take into 
account both the influence of rotation and magnetic fileds. 
The second string in ^ is responsible for the heating and 
cooling of the matter in CZ due to generation and dissipation 
of large-scale magnetic fileds and differential rotation there. 

The angular momentum transport in convective shell is 
described with, 
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The turbulent stresses, Ti 



(r,e)<t>, 



describe the effective redis- 



tribution of the angular momentum in rotating convective 
shell due to convection. The mesure of rotational influence 
on convection is the Coriolis number, Q.* = 2r2oT"c; here 
f2ois the basic rotation rate of a star, and TcIs the typi- 
cal convective turnover time. The parameter /3 mesures the 
influence of magnetic filed on convection. It is defined as 
(3 = \B\ I (v^ivrpuc), where Uc is rms convective velocity. 
Both Tc and are taken from the reference solar interior 
model by Stix(2002). Functions, <E>(^j|), ^{1,2) are respon- 
sible for the influence of rotation on the turbulent viscos- 
ity and on the A effect repectively. They are given in pa- 
pers([5) |9|). The effect of magnetic fileds on viscosity and 
on the A effect (which is for breivity call as A-quenching) 
is approximated with help of 'I'i,2 and with ip(v,H) respec- 
tively. They are defined as follows. Its definitions can be 
found in appendix. In case of the slow rotation we have 



*{1,2) («*,/?) = ^(1.2) (P), Av,H) (O.' 

^{1,2) are given in ... and i^(i.2)- in ., 



P) = K (/3)(j 2)"^here 
The Cv is a paramer 



to adjust the differential rotation law in shell in more close 
argreement with helioseismology results. We use cv ~ 1.3 
to get nearly radial profiles of angular velosity isolines in 
the main part of CZ (see below). — — is the effective 

XT 

Prandtl number and xt, i^t - typical eddy heat conductiv- 
ity and eddy viscosity respectively. In simulations preset- 
nted here we used P^ ~ 0.3,0.1,0.05. These gives examples 
for two different regimes. The effects of A-quenching pre- 
vails the Lorentz force for = 0.3 and the opposite case 
is for Pi, — 0.05. Aslo we have to introduce the parameter 
Fl < 1 to tune the power of the large-scale Lorentz force. Its 
necessity is determined by shortcomings in dynamo model. 
Particularily, in our model the relation between the typical 
strength of toroidal and poloidal components of the large- 
scale magnetic field is Br/Bp ~ 50 which is less than can 
be expected from observations, where Br/Bp > 100. 

We use the following dynamo equations in a spherical 
coordinate system, 
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More details about the given dynamo equations can be found 
in ... Below, we make some commentaries about (3j4|. The 



power of generation sources due to alpha effect is character- 



ized via Cq 



770 



where 770 and ao are the maxima of 



choRq 

radial distributions of magnetic diffusivity and alpha effect 

in CZ; ajv/ ~ -77^ is a mixing-length paramer, which relates 
Hp 

the pressure scale hieght HpSnd typical mixing length of 
convective flows £c; Pm relates the eddy diffusivity and eddy 
heat conductivity tjt = PmXT- To match the period cycle 
we used = 0.05 in what follows. Functions 4>n,\\, .4nand 
'^cm surve to determine the dependence of the turbulent 
drift, dissipation and generation of magnetic fields on the 
Coriolis number. They are discussed in ...[]. The magnetic 
quenching function for magnetic diffusivity and transport 
effects is ■)/',,. The a quenching is described with (j)a. Both 
are given in appendix. The eddy viscosity is determined with 
help of mixing-length relations for rms convective velocity 

Ucic 
U-c, XT = ^5—1 



Cv dr ' 



2.1 Boundary conditions 



The system ( 2|3|4|1 1 is integrate on the volume domain 
which is confined on the radius in r = [0.715 — 0.96] Rq 
and on the polar angle - from pole to pole. For the an- 
gular momentum balace problem we define the stress-free 
conditions both at the bottom and at the top of CZ via 

Trcft H 7, 7: -iTx = 0. For dynamo equations, at the bot- 

47rr^ sm aO 

tom of CZ we put a superconductivity condition. At the 
surface we assume the condition of the "partial reflection" 
of magnetic energy suggested previously by Kitchatinov et 
al. .., i.e., for vector potential A the vaccum conditions is 
applied and the asimuthal filed B is related with poloidal 



component of electromotive force via £0 = 



, where 



we takeS, 



■0-D(O) 



(0) 



200Gs. For the heat transfer problem, at the 



bottom we take the vanishing convective flux and at the 
surface the thermal flux is approximated an effective black- 



body radiation: Fr 



1 + 4: 



TJS 



where T, 



is the the effectictive temperature of the photosphere and 
Te is the temperature at the outer boundary of the inte- 
gration domain. The reference stratification inside domain 
is assumed nearly adiabatic. It is fixed if to specify the 
temperature Teand density peat the outer boundary. We 
use the parameters from the solar interior model by Stix, 
Te = 1.93 X lOpKand pe = 3.6 x 10"^g~cm"^ 



3 RESULTS AND DISCUSSION. 

Before discussion results it remains to fix some of param- 
eteres introduced above. The mixing-length parameter is 
specified to = 1.3. The low magnetic Prandtl number 
Pm = 0.05 helps to match the period of magnetic cycle in 
the model. Note, the maximum of radial heat conductivity 



in our model is xt ~ 2.1 x 10^^cm'^s~^so the correspon- 
dent magnetic diffusivity is about lO^^cm'^s"^. The threesh- 
old Cawas obtained via numerical trials. As an intial con- 
dition we choose the weak magnetic field with zero parity 
index (see the definition below). The dipole parity activates 
first for > 0.66. For the Prandtl number = 1 we 
did not find the long cycle within range 0.66 < C'c, < 1. 
Note, the choice Ca > 1 would mean that we are going 
to inject more energy in alpha-effect than it is available 
from convective flows. The long cycles activate after increase 
the hydrodynamical response time on the angular momen- 
tum balance pertubations by magnetic flelds. For analyzing 
purpose it is useful to introduce the following global quan- 
tities of solution. The overall magnetic activity is quanti- 
fied by integral of modulus of magnetic flux contained in 
CZ, Fm = Jo Jr" \B {r,e)\r sin edrde The parity index - 

P = -prrjT Fu^' where E^^\ iS'^' the integrated over 

-j- F^ ' 

convection zone energies of the quadrupole and dipole com- 
ponents of magnetic flelds, respectively. 

The main subject of the paper concerns with the influ- 
ence of the shear contributions to alpha effect on the long- 
term variations in solar dynamo. The system ( 2|3|4|1 \ forms 
the solid self-consitent basis for investigation this problem 
with numerical simulations. Previous studies show that in- 
teraction of the LSMF and differential rotation is, perhaps, 
the most important factor that is responsible for excitation 
of the secular magnetic activity cycles in the non-linear stel- 
lar dynamo. Coupling the LSMF generation and angular 
momentum transport is described by 1 2|3|4 1. The changes 
of the heat transport in CZ due to LSMF activity and differ- 
ential rotation result to modifications both in the dynamo 
processes and in the angular momentum balance. The in- 
flunce of these modifications to the secular magnetic activity 
variations needs the separate study. We will not address this 
problem here. At the end we show some preliminary results 
about influence of the LSMF and differential rotation vari- 
ations on the changes in solar irradiance and radius. Just 
below we consider results where the heat transport problem 
([T| is fixed at the beginning of time-stepping procedure by 
solving without regards for the LSMF. 

The Fig[l] shows some results of numerical simulations 
in the case of Ca = 0.85. On the left side of the Fig[T]we 
put the results concerning the variations of the integral mag- 
netic flux FAf(top), parity P(middle) and latitudinal shear 
(bottom). The shear is quantifled by the difference of an- 
gular velocity between equator and poles. This difference is 
averaged over N and S-hemispheres. There with solid line 
we draw the result of simulated evolutions with regards for 
shear-induced alpha effect. The dashed line shows the same 
but with the discarded contribution of shear to alpha ef- 
fect. In the picture we flltered out the short cycle to see 
the long-term behaviour more clearly. In both case we see 
the long-term variation of the magnetic activity. They are 
caused (similar to ( |19[ )) primerely by the global interaction 
of LSMF and differential rotation. It is clearly seen that pat- 
tern of the long-term variations is much more complicated 
if the shear-induced alpha effect is allowed for. Though this 
effect is certainly not the principal factor in exciting the 
long-term variation. In particular, for this model the most 
important factors driving the long-term variations are the A- 
quenching, the large-scale Lorentz force and the non-linear 
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Figure 1. Variations of magnetic flux and parity (lelt) and the Maunder diagramms for toroidal magnetic fields and angular velocity 
variations (right). 



magnetic diffusion. The long-term variations disappears if 
one of them is excluded. 

On the right side of the Figjl] we draw the butterfly 
diagramm of toroidal LSMF near the surface for time inter- 
val about 200 years (top). The amplitude of magnetic field 
changes is 400 Gs. The corresponding diagramm of the tor- 
sional variations of angular velocity at the surface for this 
period is given at the middle. The variations is measured by 
deviations of the current value of angular velocity from the 
moving mean quantity obtained with averaging on the pe- 
riod of cycle. The frequency of torsional variations is twice 
to the basic frequency of magnetic cycle. The magnitude of 
the deviations of roatational velocity is about Am/s. The 
diagramm of the long-term variations of angular velocity 
spanning the whole interval of simulations is at the bottom. 
The short-term torsional oscillations are filtered out there. 
The negative deviations of velocity are shown with dashed 
contours and positive ones are shown with solid lines. The 
amplitude of seqular modulation of surface angular veloc- 
ity is about 40 m/s. According to the given diagramm the 
relatively accelerated (or slow down) zones occupy the high 
latitude of the solar surface and they drift to poles in course 
of grand activity cycle. On the picture some isolines span 



from pole to pole what can be interpreted as the angular 
momentum exchanges between hemispheres. This confirms 
the global nature of the long-term variations demonstrated 
by our model. 

At the end we discuss briefly the possible impact of 
the long-term variations of magnetic activity and differen- 
tial rotation to the heat transport and hydrostatic balance 
in SCZ. The heat transport is described with eq.([T]). The de- 
scription of hydrostatic balance with regards for LSMF and 
differential rotation is given in , .The contribution of mag- 
netic field to irradiance output is estimated by the outflux of 
the LSMF's energy from the dynamo region. This gives the 
maximum amplitude of variations in solar irradiance if we 
assume that the energy of escaping field is completely con- 
verted to radiation. The radius variations were estimated 
from the changes of the zero mode of the gravitational po- 
tential. On the Fig|2] we show of the long-term variations 
of the integral parameters describing the integral magnetic 
fiux and parity (top), the latidudinal shear and the relative 
radius variations (middle), the relative variations of luminoc- 
ity and the heat outflux (bottom). The short cycles were 
fltered out there. The panel in the middle of Fig|2] clearly 
shows that variations of the latitudinal shear and solar ra- 
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Figure 2. Long-term variations of the integral magnetic flux (top, solid line, left scale), parity (top, dashed line, right scale), latitudinal 
shear (middle, solid line, left scale), relative radius variations (middle, dashed line, right scale), relative variations of luminocity (bottom, 
solid line, left scale) and the relative heat outfiux variations (bottom, dashed line, right scale). 



dius go in phase, while they vary slightly ahead the changes 
of the magnetic activity. In our model the delay is about 
100 years. It worth to note that in the short, 11-th year cy- 
cles the radius varies nearly in phase with magnetic cycle 
with relative amplitude about W~^Rq. The amplitude of 
secular variations of radius is larger but it still rather small, 
1.5 X IQ-'^Rq. This model seems to confirm our previus re- 
sult. In .. we concluded that the primary source of the solar 
radius changes could be the modulation of the centrifugal 
forces by magnetic activity (via modulation of differential 
rotation). In according to simulations the secular modula- 
tion of luminocity due to the escape of the magnetic energy 
from the dynamo region is about 2 x 10~^£q. This is the 
maximum possible estimate for the LSMF input to radiation. 
The seqular variations of the heat outflux in the model are 
rather small, 10^*^0. They are in anti-phase relation with 
magnetic activity. 

In conclusion we would like to summaries the results of 
the paper as follows. The differential rotation modifies the 
standard alpha effect changing its amplitude and inducing 
the additional mean electromotive force that is perpendecu- 
lar to original magnetic field. The dependence of the alpha 
effect on differential can complicate the long-term variations 
in non- linear dynamo. The main shortcomming of the model 
presented is that we ignored the helicity conservation law. 
Some preliminary results given in .. show that the helicity 
conservation could controll the largest time scale of long- 
term variation in the non-linear dynamo. The present model 
demonstrates the secular variations of radius. The secular 
maximum of the radius corresponds to epoch the maximal 



differential rotation and it foregoes the secular maximum of 
magnetic activity. The amplitude of of secular variations of 
radius is 1.5 x 10~^i?Q. The corresponded variations of solar 
irradiance are about 2 x 10~^>Cq . Given results suggest the 
futher work in direction We hope that the futher develop- 
ment of the dynamo theory and the model elaboration will 
allow to make more solid conclusions about influence of the 
nonlinear dynamo pr changes of the global parameters of 
the Sun refinement 
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